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We explore the physical properties of a unified microscopic theory for the coexistence of supercon- 
ductivity and charge density waves in two-dimensional transition metal dichalcogenides. In the case 
of particle- hole symmetry the elementary particles are Dirac fermions at the nodes of the charge den- 
sity wave gap. When particle-hole symmetry is broken electron (hole) pockets are formed around the 
Fermi surface. The superconducting ground state emerges from the pairing of nodal quasi-particles 
mediated by acoustic phonons via a piezoelectric coupling. We calculate several properties in the 
s-wave superconducting phase, including specific heat, ultra-sound absorption, nuclear magnetic 
relaxation, thermal, and optical conductivities. In the case with particle-hole symmetry, the specific 
heat jump at the transition deviates strongly from ordinary superconductors. The NMR response 
shows an anomalous anisotropy due to the broken time-reversal symmetry of the superconducting 
gap, induced by the triple CDW state. The loss of lattice inversion symmetry in the charge density 
wave phase leads to anomalous coherence factors in the optical conductivity and to the appearance 
of an absorption edge at the optical gap energy. Furthermore, optical and thermal conductivities 
display anomalous peaks in the infrared when particle-hole symmetry is broken. 
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INTRODUCTION 



The quasi two-dimensional (2D) transition metal 
dichalcogenides (TMD) 2H-TaSe 2 , 2H-TaS 2 and 2H- 
NbSe2 are layered compounds where s-wave supercon- 
ductivity coexists with a charge density wave (CDW) 1 ^ 
at low temperatures, and whose transport properties 
are highly anisotropic in the high temperature CDW 
phased. There is a vast literature reporting anomalous 
effects in the CDW phase, including, non-linear Hall 
effect, anomalous impurity effects in the superconduct- 
ing (SC) phase 4 ", stripe phases^, and different regimes 
of commensurability^. Recent angle resolved photoemis- 
sion experiments (ARPES) reveal that the quasi-particles 
of TaSe2 have a marginal Fermi liquid (MFL) lifetime^. 
This scenario becomes more exciting by the verification 
that some of the physical properties of TMD, such as the 
linear growth of the normal resistivity with temperature 3 , 
and the strong anisotropy in the in-plane and out-of- 
plane transport are similar to the same properties in 
the high temperature superconductors (HTc). HTc do 
not show a CDW gap but a d— wave pseudo-gap coex- 
isting with the superconducting phase. In both cases, 
the transport and thermodynamic properties are weakly 
dependent on the application of external fields in the 
normal/pseudo-gap phase, and strongly dependent on 
them in the superconducting phased. Furthermore, the 
application of pressure in TMD favors the superconduc- 
tivity and suppresses the CDW phased, in a close anal- 
ogy with the HTc phase diagram of temperature versus 
doping. Differently from the HTc, however, the TMD are 
very clean materials. The anomalous TMD properties are 
sample independent and can help to clarify the physics 
behind a whole class of low-dimensional superconductors. 



The interpretation of the experimental data in TMD is 
however still very controversial. Within the Peierls the- 
ory, the CDW gap formation in ID systems is usually 
due to nested Fermi surfaces. In 2D systems, the nesting 
is not perfect and some parts of the Fermi surface may 
not be gaped. Early band structure calculations^ indi- 
cated that the T centered sheets (Si) are nested with the 
K centered ones (Sn) by the Qi (i — 1,2,3) wavevec- 
tors of the triple-CDW (see Fig. 1). The value of 
the CDW wave-vector, |Qj| ~ |rif, measured by neu- 
tron diffractioniiii^, and recent scanning tunneling mi- 
croscopy (STM) experiments^i^i^ confirm the plausi- 
bility of a nesting scenario. An alternative theory pro- 
posed by Rice and Scott 16 is based on a Fermi surface 
independent CDW mechanism, where the CDW wavevec- 
tors connect the saddle points (indicated in Fig. 1 around 
\TK) of the transition metal d— bands, generating a log- 
arithmic divergence in the electronic susceptibility. How- 
ever, the saddle point energy in NbSe2 is too large (~ 50 
meV) in comparison to the CDW ordering thermal en- 
ergy ksTcDW ~ 3 meV to allow a saddle point driven 
instabilitjiii. In TaSe2, however, ARPES has observed 
an extended saddle band along TK. This band is nearly 
flat and closer to the Fermi energy than the band calcu- 
lations predicte d 1 ^ 1 !? . As the saddle points are not well 
defined in this case, it is questionable to justify the CDW 
wave- vector measured with neutrons by some mechanism 
related to special parts of the saddle bands. More exper- 
imental studies are required to elucidate this point. 

If on one hand these arguments seem to rule out 
at least a conventional saddle point mechanism, con- 
sensus on the origin of the CDW instability has not 
been reached. STM scans at 4.2 K in TaSe 2 , TaS 2 
and NbSe2 show that the amplitude of the CDW gap 
is A C dw ~ 80, 50, and 34 meV, respectively^. The 



ability of ARPES to measure the superconducting gap 
A s ~ 1 meV <C Acdw in NbSe2, combined with the 
complete failure of ARPES to detect traces of the CDW 
gap in the Brillouin zone of TaSe2 and NbSe ^ 7 i 21 were in- 
terpreted as an evidence that the Fermi surface is weakly 
covered by the CDW. We observe that the ARPES results 
seem to be in contradiction with the STM data, and can- 
not explain the non-Fermi liquid transport in the TaSe2 
crystal. One possibility is that the ARPES data are ob- 
scured by the strong dependence of CDW gap with the 
directions of the Brillouin zone combined with the forma- 
tion of pockets in the points of the Fermi surface where 
A c (k) = (max[A c (k)] = Acdw)- Another possibil- 
ity is that the ARPES electronic dipole matrix elements 
vanish for certain states in the CDW phase due to the 
broken spacial inversion symmetry (detected in neutron 
scattering^) forbidding the observation of some bands. 

The strong resemblance of the normal CDW phase re- 
sistivity of TaSe2 with the HTo 2 ^ and the anomalous 
quasi-particle life-time, given by the inverse of the imag- 
inary part of the electronic self-energy^ ImE(fcj?,a;) oc 
t _1 +6|oj|, indicates that a marginal Fermi liquid theory^ 
should be developed as the basis of a minimal model uni- 
fying the CDW and superconducting phases in TaSe2- 
The experimental verification that ksTcDW <C Aqdw 
for all the TMD (in TaSe2 for example, ksTcDW ~ 
120 K ~ 12 meV) gives a good indication that a strong 
coupling CDW theory is required. 

One of us (A.H.C.N.)2i has recently proposed a uni- 
fied picture for the CDW and SC phases where the el- 
ementary particles are Dirac fermions that are created 
in the region where the CDW gap vanishes, leading to 
the generation of a nodal liquid. According to neutron 
diffraction studies, the inversion center of the crystal is 
lost in the CDW phased, allowing for the possibility of 
piezoelectric effects. In a system with nodal quasipar- 
ticles, the piezoelectric coupling is a marginal coupling 
from the renormalization group (RG) point of view, while 
the usual electron-phonon coupling is irrelevant under 
the RG2&. Based on a tight-binding description of the 
electronic orbital^, and on the assumption of imperfect 
nestin g be tween different Fermi surface sheets, the model 
of ref. [2J| proposes a /-wave symmetry CDW gap, with 
lobes along the saddle point directions, and six nodes 
at the points where the gap is zero (see Fig. 1). The 
proposed CDW gap is odd in the Brillouin zone, due 
to the symmetry of the electron-phonon coupling 2 ^, and 
due to the absence of inversion symmetry in the CDW 
phase, changing sign in each node. The superconduc- 
tivity emerges from Cooper pairing between the Dirac 
fermions mediated by acoustic phonons via a piezoelec- 
tric coupling. We propose that the Fermi surface is fully 
gaped by the superposition of the CDW and the s— wave 
superconducting (SC) order parameters. This model is 
able to correctly explain some of the anomalous proper- 
ties of the TMD like the marginal quasi-particles life-time 
in TaSe2 , the dependence of the normal-superconducting 
phase transition with the lattice parameters, and the 
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Figure 1: Schematic representation of the TMD Fermi surface. 
The T centered sheet (Si) is nested with the K centered ones 
{Sii) by the CDW wavevectors Qi. A CDW gap develops in 
the two sheets, except in the nodal points, indicated by the 
black filled circles. The empty circles are the saddle points. 
The thick solid line around the Si sheet is the proposed CDW 
gap. The dashed lines indicate the nodes connected by Qi. 



metallic behavior of the resistivity in the CDW phased. 

The geometry of the proposed CDW gap is similar to 
the Brillouin zone of graphite, where the nodes represent 
the points where the conduction and valence tt— bands 
cross each othe»2I. In contrast to graphite, the lattice 
inversion symmetry is broken in the distorted phase, and 
the piezoelectricity can arise. As it is usually observed 
in insulators, since metals screen the polarization fields, 
one may ask: is actually possible to find piezoelectric- 
ity in a superconductor? To answer this question, we 
should consider first that in a nodal liquid the density of 
states (DOS) goes to zero in the nodes, and therefore the 
electrons cannot effectively screen electric fields. Hence, 
one can conciliate a metallic theory (with gapless quasi- 
particle excitations) with piezoelectricity. The rigorous 
vanishing of the DOS at the Fermi surface, however, is 
not essential for the piezoelectricity to appear. It is suffi- 
cient to consider that the electrons of low lying momen- 
tum (for example, in a small pocket around the nodes) 
are "slow" enough to couple with the acoustic phonons 
of the polarized lattice. 

If the piezoelectricity and the metallic character are 
not mutually excluding, it remains the question of how 
the polarization vector affects the phase coherence of 
the condensate. The answer to this question can be 
found in the collective modes. The electromagnetic 
gauge invariance of the superconducting state is pro- 
vided by the longitudinal response of the collective ex- 
citations, that screen the electrons through a cloud of 
virtual plasmons 2 ^. Only the plasmons respond to the 
longitudinal fields and give rise to screening. Since the 
piezoelectricity involves electric fields only, it does not 
affect the phase coherence of the electrons. In a previous 



work^iL, we have shown by means of a semi-classical cal- 
culation that piezoelectricity is not only consistent with 
the stability of the condensate as it is possibly behind the 
quantum critical points (QCP) observed experimentally 
in the TaSe2 phase diagram, separating the T = com- 
mensurate phases from the stripe phase as a function of 
the applied pressure. 

The organization of the paper is as follows: in sec. II 
we introduce the model Hamiltonian of the CDW and su- 
perconducting phase; in sec. Ill we derive the supercon- 
ducting gap equation; sec. IV is devoted to the thermo- 
dynamics of the superconducting phase while in sec. V 
we calculate the acoustic attenuation rate and the NMR 
response; in sec. VI we calculate the optical and thermal 
conductivities; in sec. VII we discuss the Meissner effect; 
sec. VIII contains our conclusions. 
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Figure 2: Nesting condition in the two sheets Si and Sn of the 
TMD Fermi surface. The momentum k outside Si is mapped 
by a CDW wavevector displacement into k + Q, inside Sn. 
As the free electron dispersion £k is odd with respect to the 
Fermi surface sheets, we have £k = — £k+Q. 



II. THE HAMILTONIAN 

The nodal system is composed of two subsystems de- 
fined by the nodes of the CDW state which are connected 
by the triple-CDW wavevectors Q; (i — 1, 2, 3). It is con- 
venient to introduce the spinors 

•«*>-Ux)-(&$)- 

where +, — indicate the two nodal spaces, where 
c k(r( Ck . cr ) are creation (annihilation) operators for elec- 
trons with momentum k and spin a =t, 4- The electronic 
Hamiltonian in the normal CDW phase is composed of 
two terms, 

Hcdw = H e + H e - C . 

H e is the Hamiltonian of the free electrons in the vicinity 
of the nodes, 

= \ E ^(k)[( £ k + ek +Q! K fc 

\i,a,b.<j,i 

+ (ek-£k+ Q! )r?3 fc ] Vw(k), (1) 

where r\ v (v = 0, 1, 2, 3) are Pauli matrices that act in the 
nodal indexes a, b = ±, and £k is the free electron disper- 
sion. In our convention, r] is the identity and v = 1,2,3 
indexes the x, y, z directions, respectively. The second 
term in the Hamiltonian, H e ^ c , is the CDW exchange 
Hamiltonian between electrons situated in two different 
nodes connected by Qi, 

H e - C = Ackcj^Ck+q, +h.c. 

i,k 

= J2 A ek<i,a( k H a W(k) • (2) 
2,k,a".a.b 



where A c k is the CDW gap, with odd-parity in the nodal 
space due to the loss of the lattice inversion symmetry. 
This term arises from the scattering of the electronic wave 
function by the CDW periodic superstructure. 

Applying the nesting condition ek + £k+Q = 0, (see 
Fig. 2) in Eq. JQ, and taking the long- wavelength, low- 
energy limit, the Hamiltonian in the CDW phase reads, 

Hcdw = E *L( k ) [ v Fk±rj 3 + v A k\\r)i] ^(k), 

k,er,i 

(3) 

where and fey are the momentum components in 
the normal and parallel directions to the Fermi sur- 
face, respectively, vf is the Fermi surface velocity and, 
f A = The CDW elementary excitations around 

the nodes are therefore fermions which follow the two- 
dimensional massless Dirac Hamiltonian, similarly to the 
two-band electronic description of graphiteSl. 

The broken lattice inversion symmetry due to the 
CDW gap allows the piezoelectricity in the crystal. We 
propose that the electron-phonon coupling is piezoelec- 
tric, giving rise to a pairing of Dirac fermions in the sin- 
glet state through the triple-CDW superstructure. In 
contrast with usual Cooper pairs, whose electrons are 
paired across the Fermi surface, these pairs are formed 
by electrons located in different nodes linked by a CDW 
wavevector Qi (see Fig. 1). The pairing approxima- 
tion consists in assuming a condensate of pairs whose 
center of mass have momentum and zero spin. This 
assumption clearly violates the time-reversal symmetry 
of the superconductor order parameter A s . According to 
Anderson^, the strong insensitivity of the BCS super- 
conductors to impurities is due to the tendency of elec- 
trons to be in the state of highest possible degeneracy in 
the condensate, implying pairing each electron with its 
symmetric in spin and momentum. In such case, the scat- 
tering channels promote transitions between two degen- 
erated states, keeping the system coherent. The absence 
of time reversal symmetry should destroy the condensate 
in the presence of a very small impurity concentration^!. 
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In the case TMD, however, the CDW scattering does not 
affect the degeneracy of the condensate as far as the Dirac 
fermions living in different nodal subspaces (indexed 
by the three CDW directions i — 1,2,3) remain decou- 
pled. For this reason, we may drop the i index from now 
on. 

After tracing the phonons, the piezoelectric pair- 
interaction has the form24, 



Hp 



-fJ 



E E % a6 % cd < T (k)^ ;( _ k) 



k,k' a.b.cd 



x^ T (k')^,i(-k') 



where g is the coupling constant. The choice of the anti- 
symmetric Pauli matrix rj2 incorporates the broken sym- 
metry of the superconducting gap. In the mean field 
approximation, the pairing Hamiltonian reads, 



Hp 



where 



EE[ A wi T ( k )%°Xi(-k)+/>.c. 

k a.b 



As-.gEE^tWifi.il-k)) 

k a.b 



9 

(4) 
(5) 



is the complex superconductor order parameter. 

So far, we have discussed the problem with particle- 
hole symmetry, that is, the chemical potential [i is ex- 
actly at the Dirac point (fi = 0). In order to include 
the situation where particle-hole symmetry is broken we 
have add to Eq. (EJ a chemical potential term 



(6) 



This term introduces an electron (/i > 0) or hole (// < 0) 
pocket around the Dirac point producing a finite density 
of states. 

In order to diagonalize the problem it is convenient to 
extend the spinorial notation to the Nambu space 



*(k) 



/ V>+,T(k) \ 

W! a (-k)/ 



(7) 



with k defined with respect to the nodes. We introduce 
a new set o Pauli matrices r M which operates in the space 
(f k, I —k). Denoting r^r/^ as the tensor product between 
the Nambu and nodal spaces, it is not difficult to see that 
the full Hamiltonian is written as 

H = ^(k) [v F kj_T r/ 3 + VAk\\T Q f]i 

k 

+A sTl ri2 - piT 3 r] Q } *(k) . (8) 



Notice that the gauge symmetry of the problem ip — > 
V>e lfl , and A s e 2ie — > A s , is broken at the mean-field level. 
With this notation, the SC order parameter is given by : 



A s = - 5 ^($t(k)r 1I)2 $(k)) 



(9) 



The diagonalization of the Hamiltonian leads to four 
branches of excitations: 



±E } 



= ± 



yj {v F k ± /J,) 



2 + A?, 



(10) 



where k = k^ + {v&/vF)k\\ is the in-plane anisotropic mo- 
mentum, with k = |k|. In the normal phase, we identify 
two branches of excitations (we assume /j, > 0) : 



ztvpk + fjj (hole- like branch) 
ztvpk — fi (particle- like branch) , 



which are related to hole and particle-like pockets around 
the CDW nodes (for \x < 0, the nomenclature is ex- 
changed). The two branches are physically equivalent 
to each other, except for a constant equal to — ^ k 2fi, 
integrated in the volume o f the Dir ac cone. The optical 
gap in the SC phase is 2y7? + Af, as one can see from 
Fig. 3. 




Figure 3: Dirac fermion dispersion in the pocket with the 
opening of the SC gap for = | A s = 1 meV. Each band has 
two pocket branches indicated by the dotted and thick solid 
lines. The thin solid lines with the vertex above (below) the 
Fermi energy E = represent the hole (particle)-like branches 
of the Dirac cone in the normal CDW phase. 



III. THE GAP EQUATION 

To calculate the gap self-consistently, we use the stan- 
dard many-body Green's function method. Since Hamil- 
tonian ijHJl has a quadratic form, H = X^k 1 ^ ^ ^ , its 
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corresponding Green function in the 4x4 space is: 



G (k, iLU„ 



dre^" r (T r [*4' t ]) 



where T r is the time-ordering operator in imaginary 
time, Lo n are the fermionic Matsubara frequencies, fi = 
l/(fcsT) is the inverse of temperature, k B is the Boltz- 
mann constant, and 



W k = V F T T]- k+ A S T 1 ?72 - fJ.T 3 T) . 



(11) 



is the dispersion tensor with vpfj ■ k = vpk±V3 + ^A^H^i- 
Exploring the anti-commutative property of the Pauli 
matrices, the Green function which is systematically used 
in our calculation is: 



G K,k) 



E{/ + 2fj,v F T 3 fJ-k 



K 2 



(12) 



where 



K 2 



v\k 2 + m 2 + A 2 S = E 2 ^ - 2v F (±u)k . 



Noting that (k)*^)) is the retarded part of the 
Green function, Gp a (k, r — > 0_), we see from Eq. © 
that the amplitude of the mean-field gap is written in 
the Nambu notation as 



2A, 



E E Tr [r^a G (w„,k) 



k uj n — — oo 

Evaluating the trace yields: 



2A, 



2ttva 



<7=±1 



A 



dk 

-fc'k.tr^ 



tanh 



(13) 



where A is a momentum cut-off associated with the lin- 
earization of the dispersion close to the CDW nodes. 
For /z = 0, the gap equation is rather simple and reads, 



A s (T )ff)M =0) 



cosh 1 [cos\i[kv&vfI3 / g c 



xe 



-TrVAVF/3/g 



(14) 



where g c — 27tva/A is the zero temperature critical cou- 
pling constant. In fact, 



A S (T = 0,g, ft = 0) = 2 7 r UA i; F .g- 1 1 



9c 
9 



(15) 



Notice that for g < g c we find A S (T = 0, g < g c , fi = 
0) = 0. Hence, the // = gap equation has a quantum 
critical point (QCP), indicating that superconductivity 
occurs only above a minimal coupling g c . This is a gen- 
eral property of the nodal liquid due to the absence of 
the background Fermi sea. In a Fermi liquid (where the 
Fermi surface is large in comparison to all the other en- 
ergy scales), the Fermi sea is unstable to the formation 



of Cooper pairs between two electrons mediated by an 
attractive potential, even for infinitesimal coupling^. In 
this case, the Pauli exclusion principle of the background 
electrons plays the role of the interaction, making the 
condensate stable even in the weak coupling limi^. The 
zero temperature gap ltl5)l equals to the energy cut-off 
a = vpK in the g — > oo limit. 



A. Zero temperature analysis 

To see how the pocket affects the QCP when g ~ g c we 
analyze the gap equation in the zero temperature limit. 
At this point we introduce a more suitable cut-off, given 
by the momenta s± that define the surfaces of constant 
energy in the Dirac cone, 

a ee v 2 F k 2 = {v F s± ± fi) 2 + A 2 S = const. (16) 

This new definition of the cut-off (basically replacing A 
by s a , with a = ±) is convenient because it simplifies the 
integration, allowing us to find simple analytical expres- 
sions for the gap. This approximation is fairly reasonable, 
since the results of the model are not to be taken literally 
when fi and A s are comparable to the energy cut-off of 
the Dirac cone, a, in which case the contribution of the 
high energy states cannot be neglected. On the other 
hand, we should be warned by the fact that this new mo- 
mentum cut-off So- does not conserve the number of states 
of the normal phase. When calculating thermodynamic 
functions, the correct cut-off is A, which correctly maps 
the volume of the Dirac cone and avoids problems such 
as loosing states in the SC phase, what would certainly 
have an effect in the condensation energy. For almost 
all the applications, the results are not seriously affected 
by the details of the cut-off if the gap, A s , is sufficiently 
small in comparison to a. 

The T — gap equation becomes 



2A S 



^ y f"dk_k_ 



2ttvfva 



A s [2a-2^A2 + M 2 



—/z In 



(17) 



We rescale all the quantities by defining x — A s /\u\ and 



h{g) = 2itvfVa 



ImI 



The T — scale invariant equation is 



F(x, h(g)) = Vl + x 2 



(18) 
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Figure 4: Scaling of the zero temperature gap equation versus 
the coupling constant h(g) oc (g^ 1 — 1 ) / 1 /J. | . 

We see in Fig. 4 that Eq. itlSll has two distinct coupling 
regimes: 

(i) the strong coupling sector g > g c , where the 
marginal physics develops, with \fi\ <C A s (0,g,fi) for 
9 > 9c, 

(ii) the weak- coupling sector g < g C7 where the energy 
scale of the pocket is large in comparison to the gap [i.e. 
|/i| > A 5 (0, g, fi)] when g/g c -> 0. 

In the later, the system flows in the direction of a Fermi 
liquid state in the weak-coupling limit (g <C g c ), while 
in the former the nodes are well defined for g > g c , 
resulting in a nodal liquid description. We notice that 
the quasi-particle properties are strongly affected by the 
coupling constant g, which separates the marginal Fermi 
liquid (MFL) sector from the "Fermi liquid" one, where 
the pocket plays the role of the Fermi surface, raising the 
density of states in the nodes. 

For convenience, we denote the zero temperature gap 
A s (0, g,/j) by Ao M from now on. In the strong coupling 
limit, (|^|/Ao M <C 1), we may write Eq. ifTTjl as 

1= 9 L-A.+4-) , 

2ttv a v f \ 2A S J 

whose solution is 

Ao, ^ ^ (l + y/l + 2/xVAg) , (19) 

where Ao = A(T = 0, g, [i = 0) is given by Eq. Ill ."it . 
In the opposite limit, Ao M /|/x| <C 1, in the weak-coupling 
sector, we see that Eq. l|T8l) can be expanded in leading 
order in x, giving 

F(x,h(g)) ^ 1 + In (f)-%) 

= 1 + lnf^) -2nv F v A 9cl ~ 9 - 

= o, 

yielding, 

A 0M 9 ^ 2| M |e ft ^)-i 

= 2\fi\e 27TVFV&i9 - 1 ~ 9 ~ 1)ltirl ~ 1 . (20) 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 

9/9c 



Figure 5: Dependence of the zero temperature gap (normal- 
ized by the cut-off a) with the coupling constant g. Solid lines: 
numeric solution of the gap equation l|18fl: dotted: strong cou- 
pling approximation ([/Lt| -C Ao M ); dashed: weak coupling one 
^> Ao M ). We have set \n\/a = 0, 0.1 e 0.3 from the bot- 
tom to the top. Notice that the QCP at fx = is suppressed 
by the pocket formation (|/i| > 0). 



Although the strong coupling approximation is rigor- 
ously valid only for g ^> g c , and the weak coupling one 
for j < j C! these two approximations are remarkably 
good in almost the entire coupling range of their respec- 
tive sectors (as shown in Fig. 5) provided that \fi\/a is 
small. However, to find sensible results, one should con- 
sider that the valid coupling range of the theory is limited 
not too far above the critical coupling g c , in order to keep 
the ratio A 0AI /a small (see Fig. 5). 



B. Finite temperatures 

Let us return to Eq. ltl3|l . After some algebraic manip- 
ulation (see the details in appendix A), the gap equation 
in the strong coupling regime assumes the form 

C08h(/3A s /2) e^" tanh(^A s /2)/(4A s ) 

= cosh{irv A v F {3/g c ) e -™±*rP/9 . (21) 

The quantity tanh(/3A s /2)/A s changes very little with 
in the whole temperature interval. In a first approxima- 
tion, we can obtain the analytical expression of the low 
temperature gap by replacing the gap inside the exponen- 
tial by its zero temperature value Ao^ . This substitution 
leads to: 

2 

A a (T,g,fi) ~ - coslT 1 [cosh(TTv A v F P/gc) 

x e -7™ A D F/ 3/g e ^ 2 /3 tanh(/3A 0fl /2)/(4A OM ) 
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valid in strong coupling for small n/a. Close to the phase 
transition, Eq. ltl3|l gives 



A S (T) 



ft 



2 ^ 



_j_ 

tic 



7C(3) 



fa 



ImI/Ao m « 1 

(22) 

|/x|/A 0(U > 1, 



where £ = (T c — T)/T c is the reduced temperature and 
C(a;) is the Zeta function. The critical temperature is also 
calculated from the gap equation, 11 .'111 , in the A s — > 
limit, giving 



2/cb In 4 



Ao 



v/aT 



In 4 



|/i|/A 0M < 1 



Ml. p a(i-gc/g)M- 



, |M|/Aop > 1 , 
(23) 

where In 7 ~ 0.577 is the Euler constant. In the particle- 
hole symmetric case (// = 0), we have T c = Ao/(fes In 4) 
and A S (T — > T c ,<7,0) = 2Aot 5 /ln2 (see appendix A for 
details). 

We see that the existence of a pocket suppresses the 
QCP (T = 0) separating the normal and SC phases (see 
Fig. 5). This effect is due to the establishment of the 
background Fermi sea, which stabilizes the Cooper pairs 
for an arbitrarily small coupling. The thermal effect on 
the gap recovers the parametric phase transition with 
the coupling constant g, as displayed in Fig. 6 (top) by 
noting the presence of a minimal coupling (say, go(T, /z), 
with 5o(0,0) = g c ), below which A s (g < go,fJ>) = 0. 
The explanation can be found in the strong dependence 
of the critical temperature T c with g, as shown in Fig. 6 
(bottom). At a given non-zero temperature T, a minimal 
coupling is required to satisfy T c (g > g ) > T. 
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Figure 6: Top: SC gap A 3 (T, g, n) vs. the coupling constant 
g/gc- Solid lines: numeric solution of the gap equation ll.il : 
dotted lines: strong coupling solution (analytic). From left to 
right: k B T/a = 0.005, 0.1, 0.2 and \n\/a = 0.2, 0.2, 0.1, re- 
spectively. Bottom: SC gap vs. temperature. The scales are 
normalized by the energy cut-off a of the Dirac cone. Dotted 
line: (\fJ,\/a = 0.06, g/g c = 1.2); solid: {g,/a = 0, g/g c = 1.2); 
dashed: {\fJ,\/a = 0.06, g/g c = 1.1). 



IV. THERMODYNAMICS 

In this section, we calculate the thermodynamic func- 
tions starting from the partition function Z of the nodal 
fermions. The partition function is defined as usual from 
the original Hamiltonian (jHJ, written in a diagonal ba- 
sis of eigenstates indexed by k, 7 = ±1 (for the two 
particle-hole branches), a = ±1, and with eigenvalues 
Eg. = ±i?k, CT ^: 

Z = e-^ = tre-"* 

= e"^ 1 ^ I] E Kle-^K) 

k,a "£=0 

- e-^ 1 ^ J] (i + e -/")*W) , 

k,7,fj 

where is the thermodynamic potential. The Hamil- 
tonian includes the term A 2 s /g, in order to give the cor- 
rect condensation energy. The thermodynamic potential, 



ft = Qq + A 2 /g, is given by 
n (T) - -i^ln[2 + 2cosh(/3^ ffAl )] 

/■A 

= /- Y dk k In [2 + 2cosh(/3^ k , CTAI )] . 

(24) 

If fl n is the thermodynamic potential in the nor- 
mal phase, the condensation energy, O n (Q) — fio(0) = 
H 2 (0)/(87r) , is given in terms of the zero temperature 
critical field, H c , shown in Fig. 7. Analogously, the in- 
ternal energy, E = Eq + A 2 /g, is given by 

E Q (T) = ]T lE^nl 

k,7,tr 

(25) 
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Figure 7: Solid line: critical field H c dependence with tem- 
perature, in units of a/y/g^ for /i = and g/g c = 1.1; dotted: 
empirical law H C (0)[1 — T /TJ?]. The difference between the 
two curves is shown in the inset. 



Figure 8: Specific heat Cv /T vs. temperature for fi — 0, in 
units o{k%/g c . The jump occurs at UbT c = Ao/ln4. Dashed 
line: normal behavior in the absence of the SC gap. 



where = (e 7 ^ Ek CTf ' + is the Fermi-Dirac distribu- 
tion indexed by a, 7 = ±1. 

According to the standard thermodynamic relations, 
the specific heat is defined by 



C v =T— 



-0- 



(26) 



where S — (E — 0)/T = — is the electronic en- 

tropy due to the Dirac fermions. At low temperature, 
the gap is practically independent on the temperature. 
It is easy to check that the specific heat dependence with 
temperature in this limit for ^ — is given by: 



Cv 



T«T C 



1 



2-KVFV& j A 

6k B 



" A d££ 3 sech 2 (^ 



-A?, 



-0A, 



ttvfva 



leading to the expected exponential behavior of Cy with 
the gap. 

A more interesting result is related to the jump of the 
specific heat in the normal-SC phase transition. The cal- 
culation is given in appendix B for the weak and strong 
coupling regimes. It results in two well defined limits: 
the marginal one (|/i|/3 c "C 1), 



AC7 V 



Cn 



To 



21n4 
9C(3) 



In 4 



> 0.35, 



where the equality holds for /1 = 0; and the Fermi liquid 
limit 



ACV 



C n , 



1 



2tt 2 7C(3) 

8tt 2 



< 1.43. 



+ 



(27) 



which recovers the BCS result for !3 c \lA ^> 1. 

The jump observed in NbSe^V 5 ^ (AC/C„ - 2) is 
a good indication in favor of a conventional Fermi liq- 
uid and BCS behavior. In TaSe2, however, where the 



transport is marginal and the quasi-particles are not well 
defined in the Landau sensed! (rw < 1, where t~ x is the 
scattering rate), the picture can be very different. In the 
nodal liquid case, the specific heat jump is strongly at- 
tenuated due to the vanishing density of states (DOS) 
in the Fermi surface, resulting in the universal constant 
ACv/C n — 0.35. The plot of the specific heat displayed 
in Fig. 8 shows that the temperature dependence of the 
normal CDW phase is quadratic. As the DOS is raised by 
a pocket around the nodes, the jump grows in direction 
to the BCS value of 1.43, which corresponds to the weak 
coupling limit. However, we notice that the nodes cease 
to be well defined in the presence of large pockets. In 
this case, the pairing ansatz adopted in sec. II and the 
role of piezoelectricity in the electron-phonon coupling 
are questionable. 



V. COHERENCE FACTORS 

In this section we calculate two basic properties of the 
superconductor: the acoustic attenuation and the nuclear 
spin relaxation rate in the absence of impurities. 

A. Acoustic attenuation 

The ultra-sound attenuation results from the resonant 
absorption of the longitudinal phonons in the solid^. 
The absorption rate is proportional to the imaginary part 
of the charge susceptibilitj«2£: 



*(q) = -^i lim n 



-Im X c (q,w) 



(28) 



in the q — > limit, since the phonon wavelength is much 
larger than the typical electronic wavelength. This prop- 
erty is connected to the superconductor coherence fac- 
tors, which basically define the probability amplitude of 
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quasi-particle transitions between two states represented 
by the pairs space (k f,— k J.)SM&. These factors con- 
serve the time reversal symmetries of the interaction in- 
volved in the transition. They are usually divided into 
type I, for interactions which preserve the time-reversal 
symmetry (like in the electron-phonon coupling) and type 
II when this symmetry is broken, like in the spin exchange 
interaction. The charge susceptibility is defined in terms 
of the time ordered charge density correlation function. 
All the correlation functions used in this article are de- 
fined in appendix C. Using the spinor defined in Q, the 
charge density operator is given by: 

p(q) = ^4 CT (k-q/2)Vw(k + q/2) 

k.tr.a 

= * f ( k - q/ 2 ) nVo9(k + q/2) . (29) 




Figure 9: Temperature dependence of the acoustic attenua- 
tion rate normalized by the normal phase rate. Solid: our 
model (jit = and g/g c = 1.1); dashed: BCS model. 



We define: G+=G (k + q/2,wj„ + iu) and G-=G 
(k — q/2, iuj n ), so that the electronic charge susceptibility 
reads 



1 



\ " 



X c (q,iw) = -Tr ^ 



G+ T3770 G- r 3 r] . 



(30) 



It is convenient to define the gapless Dirac fermions 
dispersion by tk = v F k, and the quantity e Q = 



v 2 T (k 2 + q 2 ). After evaluating the trace and the sum 
over the fermionic Matsubara frequencies, the imaginary 
part of the susceptibility reads: 



ImX c (q,w^ 0) = 



CO 



V F 



£q/2 ?™A 



X ■ 



A 



(Tfl 



dkk 



E 

a=±l 



q/2 



dn{E 0t<JIJ ) 
(31) 



where E 0y(7 ^ = (e D + o^i) 2 + A 2 , and n is the Fermi- 
Dirac distribution. Replacing Eq. ll3T|l into Eq. l(28jl . we 
obtain the ultra-sound attenuation rate: 



1 



A? 



E 



£q/2 WAVF a £^ 1 Jo 



dee , . ^dn(E afJi ) 



(e + cr/z)- 



dE a 



The temperature dependence of a s is displayed in Fig. 9 
and shows a power-law behavior near the phase tran- 
sition. This result is compared with the BCS curve 
a./a n = 2/(e^ + 1)31. 



B. NMR relaxation 

The NMR relaxation has its origin on the hyperfine 
interaction between the nuclear spins and the electrons. 
The relaxation rate measures the nuclear spin time- 
variation along an arbitrary direction of the spin space, 
say b. The condensate exhibits no paramagnetism in 
the singlet channel, where the total spin of the pairs is 



zero. Since the Zeeman and hyperfine energies are usu- 
ally small in comparison to the gap, the only processes 
that contribute to the spin relaxation are thermally ex- 
cited quasi-particles. The inverse of the spin relaxation 
is proportional to the local magnetic susceptibility pro- 
jected along b, 



-Im^(q,w) 



(32) 



where X?(w) is given in terms of the normal directions of 



the spin space by: 



-W) X %{lo) .with 

t,j = 1,2,3 representing the x, y, z directions, respec- 
tively (see appendix C). 

Before defining the spin density operator, we must in- 
troduce the spin degrees of freedom in the spinor repre- 
sentation, Eq. 0. This is naturally done in the Balian- 
Werthamer (BW) spaced, 



*»(k) 



Ca(k) 

-io-aCj(-k) 



' <Mk) 1 
-V>L(-k) 



(33) 



V ^ T (- k ) / 

which contains an additional spin subspace 

inside the regular Nambu space, (t k, J. — k). We have 
defined a new set of Pauli matrices cr M = (<to,<?) which 
operate in this new space. The general spin density op- 
erator is 

£(q) = \ E Vi CT (k-q/2)ar'Vw(k + q/2) 

kcrcr'a 

= \Y. * t (k-q/2)a i r ?7o*(k + q/2), (34) 

keiB.Z. 

where i = 1,2,3 are the spin directions, k is summed in 
half Brillouin zone, and a, <r'=TI are the spin indexes. It 
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is not difficult to check that the Hamiltonian |JHJ , written 
in the BW space, is given by (see appendix D) 

H = ^ * f ( k ) [vFCToT rj-]i- A s (73Tl772 
kSjB.Z. 

-^croT 3 r?o]*(k). (35) 

The matrix inside the parenthesis defines the new disper- 
sion tensor lu^ for the Green function ltl2l) . G = [iw n — u) 
) . Notice that the BW Green function is very similar to 
the previous one, except for the size of the Hamiltonian 
space, which now is 8x8. 

The pairing term brings something new, because of 
the broken time-reversal symmetry of the SC phase, ex- 
pressed by the anti-symmetric property of the Pauli ma- 
trix 772 under the transposition: 77^ 6 — > rj h a a . We will soon 
explore the physical consequences of this broken symme- 
try. From Eq. IIC2J1 the spin susceptibility tensor is given 

by 



kG^B.Z. 



2J G+ <JiT Q r) G- <J 3 



tqVo ■ 



(36) 

Notice that the product OiT^riQ G ^i^~oVo —G for i = 3. 
For i = 1,2, the anti-commutative matrices i]i lead to 
a sign change in the gap term of uj inside the Green 
function, implying A s — » — A s . Thus, the i = 1, 2 (i.e. 
x, y) directions have the same coherence factors of the 
charge susceptibility, 



xL(q^) = xl 



,(q,w) = ^x c (q,^) 



This property is better illustrated in the /i 
where 



^E 



a 2 ) -w„k + w) 



[ul+El] [K+c) 2 + £ 2 ] 



(37) 
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Figure 10: Temperature dependence of the NMR relaxation 
rate normalized to the normal phase relaxation for g/g c = 1.1. 
Dashed (/J. = 0) and dotted {\fJ,\/a = 0.05) lines: NMR re- 
sponse along the in-plane directions [</?=§ case ofEq. gflt ]; 
solid (/i = 0) and dot-dashed (|/x|/a = 0.05) lines: NMR re- 
sponse along the normal c axis (</? = 0). The pocket produces 
a small Hebel-Slichter peak, indicated by the dot-dashed line. 



symmetry broken in the crystal. The calculation of the 
imaginary part of the xtz susceptibility reads: 



Uj yp I 

ImxL(q> w -> 0) = / dkk 

4:7t va Jo 



dn{E) 2A 2 e 2 /2 - M 2 



dE i k E el-fi 2 

(?k - N) ek + (|mI - ek) |a»1 



xRe 



(</ 2 -M 2 )(? k -< /2 ) 
dn(E , all ) E QtCrfl 



q/2 



e q /2e 



(38) 



where 



e k = Re 



fc q/2 



with £k = and the indexes ± representing the mo- 
mentum, ± — > k±q/2. Notice the sign difference in front 
of A 2 between the zz and the other two components. 
This difference gives rise to an axial anisotropy in the z 
direction of the spin space. This anisotropy is a conse- 
quence of the broken time-reversal symmetry induced by 
the finite momentum of the pairs, Qi, which defines the 
CDW wave vectors. This broken symmetry is reflected in 
the appearance of a spin structure oriented in the z direc- 
tion (indicated by the 03 matrix) in the pairing term of 
the BW Hamiltonian l(35|l . Therefore, we conclude that 
the z direction in the spin space corresponds to the CDW 
direction Qi in the fc-space, since it is the only rotational 



i?k = y / I^T+A 2 and e follows the definition of the pre- 
vious subsection. The Xxx,yy can be obtained from the 
substitution of Eq. (EIJ into Eq. (jHZJ|. Noting that 
X% v = X% z = Xyz — 0, the NMR relaxation rate along a 
given direction b gives 



Ti 



(b) = a 2 



— q hm 

2ir w^o 



E 



(6 2 -l) Imxt 



(39) 

In Fig. 10 we distinguish the two anisotropic principal 
directions, for in-plane b vectors and out-of-plane ones 
directed along the normal c axis. A small Hebel-Slichter 



11 



peak is formed for finite fx, but no peak is observed for 
H = 0. 

The component of the susceptibility carries coher- 
ence factors with the symmetry of the spin interactions 
(i.e. they are odd by interchanging k — > — k), while the 
xx and yy components are analogous to the charge sus- 
ceptibility [see Eq. I|37|) ]. This is easily understood by a 
qualitative argument with the aid of Eq. (|39)) . Consider 
the CDW direction Qi. The Qi direction (or equiva- 
lent^ the b = z direction for the spin, according to our 
previous discussion) affects the electronic spin correla- 
tions in the normal directions to Qi, meaning the xy 
plane. The NMR direction b = z is affected by the sus- 
ceptibility components \xx and \yy but not by the \zz 
one [see Eq. l(39jl ]. The CDW introduces an additional 
time-reversal broken symmetry to the spin correlations in 
the Qi (xy) plane, explaining why the related coherence 
factors have the same symmetry of the charge interac- 
tions. On the other hand, the planes which are normal 
to the Qi plane are affected by the Xzz component, which 
conserves the odd symmetry of the spin interactions. In 
summary, the NMR relaxation in the b = Qi direction 
(in fc-space) is therefore associated to a charge-like sym- 
metry, like in the phonon attenuation response, while the 
NMR directions which are normal to Qi have a mixed 
symmetry and exhibit a more intense response. The same 
analysis applies to the Q2,3 vectors separately. The NMR 
pattern in the fc-space results from the superposition of 
the contributions due to each vector Qi (i — 1, 2, 3) of the 
triple-CDW. As each vector Qi is rotated with respect to 
the other two by ^ and ^ (see Fig 1), if we define the 
contribution of each CDW direction to the NMR response 
along an arbitrary direction b as (b) = T-f *(£> + 6i), 
it is not difficult to verify from Eq. I|39J1 that 



1 



-3A 



2 / A d <7 



1 



2 I 7T 1 lim n ~ t Im (XL + Xlz) 



•sin Vim (xl x -X S zz) 



,(40) 



where (p is the angle that b makes with the normal direc- 
tion to the SC planes. We notice that despite the broken 
rotational symmetry of the triple-CDW state, the total 
NMR response is rotationally invariant in the planes and 
shows an anisotropic direction along the normal c-axis, 
as displayed in Fig. 10. 



VI. TRANSPORT 

In this section we calculate the optic and thermal 
conductivities of the SC phase in the clean limit. The 
transport calculation for a d-wave order parameter with 
and without d— wave superconductivity has been done by 
Yang and NayalA Here, we shall repeat the calculation 
for a CDW gap with nodes coexisting with a s-wave SC 
order parameter. We ignore the effects of scattering cen- 



ters like impurities and disorder from the CDW fluctua- 
tions motivated by the two facts: 1) the TMD are very 
clean materials, and 2) the extremely low temperatures 
where the SC phase appears in 2H-TaSe2 (T < 0.1 K), 
where conventional thermal disorder in the CDW phase 
should play no relevant role in the transport. 

The thermal current is defined by j Q = j E — , 39 
where j E is the energy current, j is the electrical current 
and /i is the chemical potential. Experimental measure- 
ments of the thermal conductivity k require zero elec- 
tric current flow in the sample, and we may assume that 
j = j E ■ The Kubo formulas for frequency dependent 
thermal conductivity, k(uj), and the optical conductivity, 
a (to), are2&: 

««M =-^f toImIlg B (q,w)+TS^(w)(7tf(w),(41) 



(TijUS) = — lira ImIIy(q,(xi) , 



where 



Sij(u) 



■— lim 

T q^O 



Imng(q,a;) 
Im Uij (q, w) 



(42) 



(43) 



is the thermoelectric conductivity (also known as ther- 
mopower) S = — AV/AT, that measures the current 
voltage AV produced by a temperature gradient AT, and 
II, IF 5 ^ and H E are respectively the electric, thermal 
and thermoelectric current correlation functions, which 
we define in appendix C. The second term in Eq. l(4*T|l 
guarantees the zero current flow condition to the charge 
carriers. 



A. Optical conductivity 

To incorporate the magnetic field into Hamiltonian (jSJ , 
we proceed with the Peierls substitution k — > k — | r 3 A. 
We assume that the vector potential A(k) is symmet- 
ric with respect to momentum inversion in the nodal 
space. For this reason, we must use the T3 Pauli ma- 
trix, which operates in the usual Nambu space. Notice 
that a given Hamiltonian density for spin \ fermions in 
the form /(k)^J (k)Vv(k) is equivalently written in 
the Nambu space as 



(>}(k) Vu(- 



/(k) 






-/(-k) 



V-t(k) 
W(-k) 



The associated matrix above is clearly T3 if / is a sym- 
metric function in k and tq if / is anti-symmetric. As the 
Dirac fermion dispersion is anti-symmetric in the cone, 
we should be especially careful with the usual Peierls sub- 
stitution, since it introduces an even term (oc T3A), which 
violates the odd symmetry of the zero field dispersion de- 
pendence with k. For a given Hamiltonian in the general 
form: 



ff = ^e k *t( k ) ToJ7 .*( k ) 
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the correct Peierls substitution involves the separation of 
symmetric (S) and anti-symmetric (A) components in k, 



e(k)r - e A (k - -r 3 A) r + e s (k - -r 3 A) r 3 , (44) 

c c 



where 
e *(k 



-r 3 A) = l\e(k- 
c 2 1 

A (k-- c r 3 A) = i[e(k- 



-r 3 A)+e(-k- -r 3 A) 

c c 



-r 3 A)-e(-k--r 3 A) 
c c 

Applying this procedure to the Hamiltonian (jHJ, it is 
easy to see that the magnetic part of the Hamiltonian is 

H B = -*t(k) \v F -A ± T Q ri 3 + VA-A^Tom] *( k ) , 
L c c J 

written in terms of To instead of t 3 , as one could naively 
expect from the straight substitution k — * k — | r 3 A. 
The current density operator j(k) = -cVa^ is given 

by 

j(k) = *t(k) [v F eTori 3 e ± + v A eTo7ne\\] *(k) . (45) 

The current-current density correlation function defined 
in appendix C is given by: 



IIx(q,iw) 
IT||(q,icj) 



P 



-Tr ^2 G+ T r} 3 G- t i]3, 

k,CJn. 

I 

-Tr ^2 G+ r ?7i G- T r?i, 



where _L and || are the normal and parallel directions to 
the Fermi surface for a given node (see Fig. 2). Ap- 
plying the Kubo formula i!42fl to the imaginary part of 
the correlation functions above, we find that the optical 
conductivity is separated into two parts: the Drude term, 




and an extra term due to the interband excitations of the 
Dirac fermions, 
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where 



4A 2 



U3* 



4/i 2 



(48) 



and 



a. 



I t-o | - fJ-)E a ^ ± (|i/ Q | -(- /i)_E ,-^| 



with ^o,±u — (l^ol ± A*) 2 + A 2 . In order to calculate the 
the parallel component cry we just have to replace vf by 
i>a- For fi = 0, the interband conductivity is given by: 



2u A 



1 



4At 



or 



1 - 2n 



-2A, 



(49) 

The conductivity ijlTj) is considerably simpler in the 
normal CDW phase. Setting the gap A s to zero, we have 



ui/2 and -Eb,±£i 



leading to: 
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(50) 



Analogously, the Drude part of the conductivity becomes 

VFe 2 fl 



DC 

'±CDW 



8v A 



de e 



x sech 2 /3 



T^O Vf& 



2v A 



S(uj) x 



21n(2) 





l/'l 



for fi = 
for ^ ^ . 



(51) 



Notice that in the absence of SC we find that <r DC (T — > 
0) is constant and proportional to fi. In the SC case, Eq. 
ijlfil) shows that a DC (T — > 0) vanishes independently of 
the pocket size (as shown in Fig. 11). The presence of a 
Drude conductivity, S(uj), results from an infinite 

electron mean free path due to the absence of scattering 
centers. If we consider that the electrons in TaSe2 have 
a finite scattering rate, T = 1/t(w)2I, the Drude peak 
will be broadened around u> = Q. The normal transport 
in the presence of an order parameter with nodes (like 
the CDW, as in our case) in the dirty limit, is given in 
Ref. 

Photon absorption involves quasiparticle excitations 
and results in the formation of in-phase currents with the 
electric field^. The absorption rate is therefore propor- 
tional to the real part of the conductivity. In conventional 
superconductors there is no absorption at T = in the 
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Figure 11: Temperature dependence of the Drude conductiv- 
ity integrated in w, for g/ g c = 1.1 and n/a = 0.1. o in units 

of «Fe 2 a/(2«A). 



infrared region where the photons with energy to < 2A S 
cannot break a Cooper pair. At finite temperature, the 
excitation channels are gradually recovered and photons 
with energy smaller than 2A S have a finite probability of 
being absorbed. We should stress that the coherence fac- 
tors in those superconductors (say, BCS type) are finite 
only in the dirty case, where the processes conserve en- 
ergy but do not conserve momentum. The first important 
distinction between traditional BCS superconductors to 
the ones discussed here is the presence of two bands, re- 
sembling the spectrum of small gap semiconductors (see 
Fig. 3). In the nodal liquid superconductor, made out of 
Cooper pairs of Dirac fermions, the absorption process 
comprehends the excitation of an electron from the lower 
to the upper band, transferring energy equal to the the 
photon energy w but with no momentum transfer. In 
the situation where the lower band is completely filled 
(fi = 0), there are no thermal channels of quasiparticle 
excitations (since the thermally excited electrons occupy 
the upper band, where there is no absorption due to mo- 
mentum conservation) and the photon is absorbed only 
when its energy is sufficient to break a pair (u> > 2A S ), 
producing quasiparticle excitations directly from the con- 
densate (pair breaking channels). When the system ex- 
hibits particle-hole symmetry, the absorption is indepen- 
dent of the temperature in the infrared for u < 2A S . 

The second important distinction between the TMD 
and BCS superconductors, is that the optical conduc- 
tiv ity show s an anomalous absorption edge in u — 
2y7? + Aj (see Fig. 12). This energy corresponds to 
the optical gap of the two bands shown in Fig. 3. The 
presence of the edge is a consequence of the broken lattice 
inversion symmetry in the CDW phase, which affects the 
coherence factors of the infrared conductivity. When the 
particle-hole symmetry is lost by shifting the chemical 
potential from the vertex of the Dirac cone, new ther- 
mal channels of quasiparticle excitations emerge, giving 
rise to an absorption peak in the infrared. To see this 
effect, we illustrate in Fig. 13 the thermal excitation 
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Figure 12: Top: optical conductivity a± vs. frequency. For 
|/x|/(2A s ) = 0.9: dashed line (T = 0) and solid (fc s T/(2A s ) = 
1.2); dotted line: \i = and T = 0. Bottom: optical conduc- 
tivity a± vs. temperature, for g/g c = 1.1 and \n\/a = 0.1. 
Dashed lines: 0.4A 0m < u> < 1.4A 0m ; solid: lu = 2.3A 0m ; 
dotted: 2.8Ao M < u < 4Ao M . In both plots, a is in units of 
ufc 2 / (2va) 



process of the hole-like branch, where the photons with 
energy smaller than 2\/j,\ are able to promote the ther- 
mally excited electrons occupying the empty states on 
the top of the lower band to the upper band. As in 
the case of superfluid He 3 , the superconductor is an elec- 
tronic liquid composed of two "fluids", where there is a 
one-to-one correspondence between the excited states in 
the SC and in the normal phases. The thermal fluc- 
tuations promote electrons from the condensate to the 
empty states above the pocket Fermi surface of the hole- 
like branch. The optical channels of absorption through 
the thermally excited electrons are therefore limited to 
the window |w| < 2|/z| (in the clean limit), as shown in 
Fig. 12 (top) and Fig. 13. 

The temperature dependence of the optical conduc- 
tivity, displayed in Fig. 12 (bottom), shows a clear 
distinction between the two absorption channels. The 
dashed lines represent the thermal channels, which van- 
ish at T — 0. The dotted lines indicate the pair breaking 
channels. These channels depend on the number of elec- 
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Figure 13: Schematic representation of the photon absorption 
process in the channel of thermal excitations of the conden- 
sate, within the absorption window \w\ < 2|/i of the hole-like 
branch (see Fig. 3). uj is the photon frequency, T represents 
the thermal excitations and /i indicates the Fermi level. 



trons in the condensate and are more effective as the 
temperature is reduced. The solid line in the same fig- 
ure represents a pair breaking channel which is abruptly 
suppressed by lowering the temperature . This is un der- 
stood by noting that the optical gap 2-y//i 2 + Af(T) [see 
Fig. 12 (top)] displaces the absorption edge towards the 
ultraviolet as the temperature is reduced. In this situ- 
ation, we expect that some of the absorption channels, 
at a given energy slightly to the right of the edge, will 
be abruptly suppressed if the temperature is sufficiently 
reduced, i.e. if the edge is sufficiently displaced to the 
right in Fig. 12. 



1. Spectral weight 



in the AC sector have exactly the same area. This be- 
havior is depicted in Fig. 14 for different values of fx, 
showing an anomalous suppression of the Meissner spec- 
tral weight at low temperatures for small fi. A superficial 
analysis would indicate that there is no spectral weight 
due to the condensate and therefore the superconductiv- 
ity is not stable. This analysis, however, is incompatible 
with the thermodynamic verification that there is a finite 
zero temperature critical field H c (0) (see Fig. 7), result- 
ing in a finite condensation energy. The origin of the 
problem has connections with the spectral weight shift 
from the high to the low energy states of the band as 
the temperature is reduced, which has been observed ex- 
perimentally in TaSe^- In this compound, part of the 
spectral weight around 60 meV (~ of the order of the 
cone cut-off) at 300 K is displaced towards the infrared at 
temperatures of the order of the SC phase transition. Ap- 
parently, the opening of the gap "attracts states" beyond 
the cone approximation. In the lowest order, the non- 
linear states in the CDW spectrum yield ek oc (k — ^A) 2 . 
These states are the only ones that contribute to the dia- 
magnetism, which results from terms oc A 2 in the energy. 
We conclude that the cone approximation excludes the 
"diamagnetic" states of the band, and for this reason the 
/-sum rule is not able to correctly incorporate the dia- 
magnetic spectral weight, specially at low temperature, 
where the contribution of the high energy states is more 
pronounced. The zero field properties which are not di- 
rectly related to the Meissner effect, however, are not so 
sensitive to the absence of the high energy states and give 
satisfactory results within the cone approximation. This 
analysis is confirmed later in sec. VII, when we discuss 
the Meissner effect in the London limit. 



According to the /—sum rule one should have 

,2 



<j{u>) dto 



nne 
2m 



(52) 



and therefore, the area "under" the curves a DC + a AC is 
conserved in the normal and in the SC phases. In the SC 
phase, however, there is a "missing" area in comparison to 
the normal phase. The difference between the two areas 
corresponds to the uj = spectral weight, responsible for 
the diamagnetic supercurrents in the Meissner effect 3 -. 
This part of the spectral weight (which properly defines 
a superconductor) depends on a different order of limits 
between w and q, and does not appear explicitly in the 
calculation. Thus, a required condition for superconduc- 
tivity is 



DC 



AC 



o 



From now on, we call the difference between the n and s 
areas as the Meissner spectral weight. 

It is not difficult to see that for \i = at zero tem- 



perature we have a 



DC 



r DC 



0, and that the curves 
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Figure 14: Meissner spectral weight A as a function of temper- 
ature. Curves drawn for < \y\/a < 0.15, from the bottom to 
the top, in fixed intervals of 0.03. A in units of i>Fe 2 a/(2i>A), 
with g/g c = 1.1. 
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B. Thermal conductivity 



The energy current is a conserved quantity denned by 
the non-diagonal components of the momentum-energy 
tensor T l , defined as^ 



dC 



(53) 



According to the usual relation H — gS^rf do^ ~ £, the 
Lagrangian associated to the Hamiltonian JHJ m the real 
space representation is 



C = ^(x) [icT 3 rj Q do - ihv F T ri 3 d 3 

-iVATariidi - A s ri?72 + /ir 3 ?7o] $>(x) , (54) 



where cdo = id T with r as the imaginary time. The 
conserved energy current j E (x) = cT l gives 



fix) 



dC 



= ^(x) [w_FT ?73e_L + WAT ?7ie||] d T ^(x) 



or equivalently 

j B (q,r) = -^*t( k _ q / 2jT )[ WFTo% e i 

k 

+v A T Q r n e ] \] w k+q/2 *(k + q/2,r) , (55) 
where the time-evolution of the Dirac fermions \& is 
*(q,T)=e- rS -*(q), 

with co defined in Eq. (JHJ. 

We are interested in the diagonal components of the 
current-current polarizations Hf E = H EE and II f x = II ;f 
given by: 



nf B (q^) = ^Tr ^ G+ TQ773 ^+ G- t 7? 3 w_ 



T0V3 



We find that 
lmnf s (0,cj) 



Vp LO 
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5(u) J2 f d 
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E <T>H 
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2 ye_ 
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2 1 
s I I 



(56) 



and 
Imllf (0,w) 



21 I ft V 

4^ Wl 1 "-^ 



[Eo,n — Eo,-n] 



X@- N^o,,,) - n(-£ ,- M )] 



xe(H-2v^+Af), (57) 



where i/o and -Eo,o> are defined as in Eq. 114 7J1 . In con- 
trast to the thermal polarization, the thermoelectric one 
does not have a Drude part. The thermal conductivity 
follows from a straightforward substitution of the pre- 
vious results H47|t. ij5fill and l|57|) into the Kubo formula 

Let us analyze these results for \i — 0. We have 



-Imnf (0,w 



VA JA 3 



E A 



X 



dn(E) 



'2va~ 



(?) 



dE 

- / 4A2 
1 + 



(w - 2A S ) , 
(58) 
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and 



-ImlF^O.w) 

to 



2v A 2 



4A 



(59) 



H-2A,). 



Replacing Eq. gjjj and |2D| into Eq. (021, the /x = 

thermopower yields 



1 ImITf(0,w) _ uj 
T Imn_L(0,cj) ~~ *MT 



(60) 



Substituting Eq. (O, ^ and (ESI into (UJ, we find 
that the only contribution comes from the Drude term 



Up 



A 2 \ dn(E) 
E 2 ) dE 



where " A( 



T , dEE 3 

vaT J As 

(61) 

k" 1 - - = for zero /i. When the system exhibits 
particle-hole symmetry, the exact cancellation of the in- 
terband contributions to the thermal conductivity is due 
to the fact that the total heat carried by a particle-hole 
pair is zero. The argument is the following^: the inter- 
band excitation process involves the annihilation of an 
electron with negative energy in the lower band, and the 
creation of a particle with positive energy — £"k + lo = 
+E k in the upper band, where to is the photon energy 
and — E'k is the energy of the annihilated electron. De- 
stroying a particle with negative energy, momentum k 
and charge e is equivalent to create a hole with momen- 
tum — k and charge — e at the energy cost +E*. The en- 
ergy current carried by the quasiparticle formed by the 
particle-hole pair is kEk + (— k)(i?k) = 0. On the other 
hand, the charge current is finite, ke + (— k)(— e) = 2ek, 
explaining why the quasiparticles are able to transport 
charge but not heat when the pocket is absent. 

When the particle-hole symmetry symmetry is lost, the 
thermal current due to the pair breaking channels is equal 
to i?k,-/j(k) + (i?k,^)(— k), or equivalently to — 2/ik in the 
normal CDW phase, when the ground state electrons are 
promoted to the upper band. As a second effect, the 
thermal channels of quasiparticle production give rise to 
an infrared peak for \u>\ < 2|/i| as shown in Fig. 15 (top), 
analogously to the optical conductivity. In contrast with 
the charge transport, however, the amount of heat carried 
by the quasiparticles is of the order of the pocket energy 
and vanishes at \i — 0. The temperature dependence of 
k is shown in Fig. 15 (bottom). The solid lines represent 
the thermal channels of quasiparticle excitation, while 
the dotted lines indicate the pair breaking channels. As 
in case of the optical conductivity, some of the latter 
channels which are slightly above the optical gap energy 
luq = 2yJ /i 2 + A 2 are suppressed at low temperatures 
(see Fig. 15 ). At T — the thermal conductivity is zero 
for \u)\ < u>o, and infinity for \lu\ >ujq. 

Let us verify the normal CDW properties (A s = 0) in 
the transport. The thermoelectric spectral function l|57JI 




Figure 15: Top: thermal conductivity xT vs. frequency. kT 
is in units of (2vf /%)Aj and lu in units of 2A S , with fi = 
2.2A S . Dashed line: T — > limit; dotted (k B T = |A S ); solid: 
(k B T = A s ); dot-dashed: (k B T = §A S ). Bellow: thermal 
conductivity dependency with temperature. We have set n 
in units of «Ffcsa/(2uA), g/g c = 1-1 and \/J./a — 0.1|. Solid 
lines: 0.4A Om < uj < 1.4Ao M ; dotted: 2A 0fl < to < 4A 0m . The 
dot-dashed line is the Drude thermal conductivity integrated 
in ui, with units of ^VFksa 2 /(2«a). 



is given by 



Imnf CDW (0, W ) 



\ ( M 




\~2" 





"Hi ^ +/! 



Comparing the expression above with the optical con- 
ductivity of the normal phase (50j) . the thermoelectric 
coefficient yields 



S±cdw(u) - i 



as in the SC particle-hole symmetric case llfiOt . Returning 
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to Eq. 115611 . and setting A s = 0, we have 



-imn 



EE 

±CDW 



(0,w) = -^- 
2va 



The thermal conductivity is given by: 

I \ DC i \ , AC i \ 

KCDW(LO) = K CDW (UJ) + K CDW (UJ) , 

where 



-DC 
l ±CDW 



M = - 



vp6(u>) 
2v A T 



E 

tr'=±l 



VpkB 

2v A p- 



■ S(ui) x 



9C(3) ,for/x = 
,for^0, 



(62) 



and 



AC i \ 
K ±CDW\. UJ ) 



,2 VF 

2v A T 



/r 
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T— >0 2 
> I" 



2v A T 
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2 



(63) 



The verification of the Wiedmann- Franz ( WF) law can 
be done in two cases. For u = 0, despite the optical 
conductivity is dominated in the low temperature region 
by the interband conductivity, 



(Ti CTH /(w,T^ 0) 



v F e~ { 

tanh' 

2v A 

+ ln(2)^<5H 



(64) 



the \w\/{kBT) <C 1 limit is dominated by the Drude part. 
Comparing the expression above with Eq. ll62l) for u = 
0, we see that the CDW phase obeys the temperature 
dependence of the WF law 



kcdw^T) _ 9C(3) 



T%Tacnw{0,T) 21n(2) 



(65) 



but with a particular numerical constant 
9C(3)/(21n(2)) w 7.8. Note that the order of the 
limits is essential, otherwise, ocdw is dominated by the 
interband term in the |w|/(fcsT) ^> 1 limit, 



lim a C Dw{u, 0) = 

w^O ZVA 



violating the WF law. We should stress, however that 
this relation is typically valid in the DC limit oj — » 0, 
which is well defined for j3\u\ <C 1 but not for j3\u\ ^> 
1. This is easily seen by noticing that at T = the 



quasiparticle excitation energy ui do not have a scale and 
the AC and DC sectors cannot be distinguished. For 
finite n, it is immediate to check that the WF relation is 
verified exactly as in a metal, 



li m k cdw(0,T) 
t^o TacDw ■ 



i0,T) 



VII. MEISSNER EFFECT 

The non-local electrodynamics is described in the Lon- 
don limit, where the vector potential function A(k) — > 
Ao = const, under the assumption that the field varia- 
tions are slow in comparison to the coherence length £. 
In this limit, the current j and the vector potential obey 
the London equation 

(ji)=j? M + Q ij A j , 

valid in the Coulomb gauge k • A = 0, where j is the 
current due to the momentum of the pair center of mass. 
For all purposes, we neglect this effect and consider only 
the response to the magnetic field. 

In order to calculate the London kernel Qij, instead 
of writing the current density operator ll4*5l) , we propose 
a more general procedure, extending the CDW band be- 
yond the cone approximation. As in sec. II, we start from 
a CDW Hamiltonian written in terms of an extended 
band 

H C dw = *^( k ) Ms + Ack»7i] *a(k) , (66) 

k,cr 

where ek and A c k are any anti-symmetric functions 
with respect to a given Fermi surface node. 

Introducing the magnetic field through the modified 
Peierls substitution J44II . the series expansion of e(k — 
f t 3 A) in powers of A is separated into symmetric and 
anti-symmetric terms in k, 

e(k--r 3 A) = [e(°)(k) + e( 2 >(k) + ...]T + 

[e«(k) 



5< 3 >(k) 



c 2 V c/ 



up to second order in A, where d l = defines the 
momentum derivatives and repeated indexes are to be 
summed. The same applies to A c (k— |t 3 A). Using the 
abbreviation k = k— -T3A, the Hamiltonian of the CDW 
+ SC phase with an external magnetic field is 

H = E* f ^) [ £ k r o?73 + A ck r ?7i 

k 

+ A s Tl772 - /iT 3 ?7o] *(k) . 



The current density operator ji (k) = —cVaH gives 
Ti (k) = *t(k) \e (ff 
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We calculate the expectation value ( j ) up to first order in 
A (see details in appendix E), and find that the London 
kernel reads, 



k a=±l 



[(ftekXfyek) + (^A ck )(^A ck )] 



xsech 



OiO,-e k + — did j A c -k 1 



k,<x/i 



x tanh 



(67) 



where = V4 + Aj and 
Ek^fj. = 



- ' A? k + ] + A 2 



is the generalized dispersion in the extended CDW band. 

The non-local properties valid in the q — > limit do not 
depend on the details of the cut-off A. For this reason, 
we are allowed to take A to infinity with no further con- 
sequences. However, the Green functions method leads 
to some spurious results in the ultraviolet if we do not 
take the full Brillouin zone into account. To see this, 
consider the illustrative case of the normal CDW band 
llfifif . After a suitable diagonalization into a particle-hole 
eigenstate basis with eigenvalues ±E£ 
we may write it into the form: 



A ck 1 



HcDW = J2 E i^( iL )^(k). 

k 

The London kernel of this problem can be derived di- 
rectly from Eq. I|67|) by setting A s = /i = 0, ignoring 
the A ck terms on it, and performing the substitution 



£k - 
has, 



ET. It is immediate to see that in this case one 



resulting in a non-zero surface term for i = j, which 
diverges in the ultraviolet for any monotonically cres- 
cent ££. The integrability of the results derived by 
this method depends on the introduction of states in 
the entire Brillouin zone. In particular, we have that 
(jf DW ) = (as expected) by assuming that the surface 
term cancels in the Brillouin zone because of its periodic- 
ity. In order to fix the spurious divergences, we follow an 
argument due to Lifshitz and Pitaevskii^. Considering 
that the kernel for A s = is zero, since no supercurrents 
are induced by the magnetic field, there is no physical re- 
sult in subtracting from the SC kernel the normal phase 
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Figure 16: London kernel dependence with temperature in 
the cone approximation (g/g c = 1.1). Plots for < \n\/a < 
0.16, from the bottom to the top, in fixed intervals of 0.02. 
0(0) — Q(A S ) in units of e 2 VFa/(2ivdvAC). 



kernel, 



&) = [Qij(A s ) - Qij(0)]Aj 



(68) 



We may consider that the kernel above correctly incor- 
porates the Brillouin zone effects, at least near the phase 
transition. 

To analyze the spectral weight behavior due to the 
Meissner effect within the cone approximation e k ~ vpk± 
and A ck ~ UA&ib we calculate the London equation in 
two limits, near the normal-SC transition and at T = 
0. Including the Brillouin zone [— 5,5] in the normal 
direction to the planes, with d the inter-plane distance, 
from Eq. © we have 



Q±(A S ) 



f3e 2 



E 

k,<r=±l 



sech 



/3i?k,o-^ 



At T = 0, the kernel gives 



.(A a )-Q±(0) 



+0 M e 2 v F 

d TTV&C 



confirming the anomalous behavior detected by the 
/—sum rule ll52ll in the optical conductivity. 

In the opposite limit, for T ~ T c , the kernel in the 
strong coupling approximation (/3 C M 1) gives 



Qx(A.) - Qj.(0) 



/3 C vf(? 
Ad it vac 



1 



M 2 f ]A» 



in agreement with the mean field result for th e penetra- 
tion depth A ± = yJc/{ATi \Q±j0) - Q±(A S )}} ocAj 1 . 

The dependence of the London kernel with /i and the 
temperature is shown in Fig 16. There is a clear suppres- 
sion of the Meissner effect in the low temperature region, 
specially when the density of states in the Fermi surface 
nodes is close to zero. As we discussed previously in sec. 
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VI, the opening of a SC gap in a nodal liquid possibly 
causes the spectral shift of high energy states beyond the 
cone cut-off a in the CDW band to the infrared. As we 
mentioned before, the spectral shift of the states bellow 
a (~ 60 meV) has actually been observed in the normal 
CDW phase of the TaSe^. More experimental studies 
are required to understand the SC phase properties in 
this crystal. 



VIII. CONCLUSIONS 

In this paper we have studied the thermodynamic and 
transport properties of a model proposed originally in 
ref. [2J| for the coexistence of a gapless CDW phase and a 
s— wave superconductor in TMD. While the lattice inver- 
sion symmetry is broken in the CDW distorted phase, as 
observed experimentally by neutron diffraction, we pro- 
pose a pairing ansatz which also violates the time-reversal 
symmetry. According to the ansatz, the pairing of the 
electrons is mediated by virtual acoustic phonons via a 
piezoelectric coupling, and the center of mass momentum 
of the pairs equals the CDW wavevectors connecting dif- 
ferent sheets of the TMD Fermi surface. This additional 
broken symmetry has dramatic consequences on the spin 
exchange interaction and produces an anisotropic NMR 
response along the normal direction to the triple-CDW 
plane. In contrast to TaSe2, the quasiparticles of NbSe2 
are well defined in the Fermi-liquid regime. The SC phase 
of the NbSe2 has been extensively studied and indicates 
that a conventional BCS description is warranted^«&2&. 

In contrast to the BCS theory, which is not critical, 
the gap equation il'M has a QCP in the critical coupling 
g = g c when the system exhibits particle-hole symmetry 
(pi = 0). When this symmetry is broken, the SC gap 
A s is strongly rescaled by \i as the coupling parameter is 
modified, and the QCP is suppressed. The scaling of the 
quantity A s //i follows two different coupling regimes: (i) 
"Fermi liquid" sector in weak coupling, for g < g c , where 
Ag/fx flows to zero as g — > 0, and (it) strong coupling 
marginal limit for g > g c , where A s /\fj,\ 3> 1. The spe- 
cific heat jump is strongly attenuated in the particle-hole 
symmetric case (where ACV/C„ = 0.35), because of the 
low density of states at the Fermi energy. As expected, 
in the Fermi liquid regime we recover the jump of the 
BCS model AC v /C n = 1.43. 

We have observed several anomalous properties in the 
transport. Unlike traditional one-band superconductors, 
the spectra for optical and thermal conductivities in the 
clean limit have an infrared peak due to the thermal chan- 
nels of quasi-particle excitation. These channels involve 
thermal intraband excitations, promoting the electrons 
in the condensate to the empty states of the pocket, at 
the top of the lower band (see Fig. 13). The absorption 
window for this channel is limited to the pocket energy 
2\fi\. A second kind of absorption channel is due to in- 
terband excitations, when a pair is broken as a result of 
the absorption of a photon. In this case, the electron is 



excited to th e upper b and, across the optical gap bar- 
rier uj — + Af . The later type depends on the 
number of electrons in the condensate and is more ef- 
fective at T — 0, except for a few channels at a given 
frequency uj a which are abruptly suppressed by the tem- 
perature reduction (say, bellow T a ) because of the optical 
gap enlargement, that is, u a < uiq(T) for T < T a . The 
thermal channels on the contrary vanish at T = with 
no exception. 

The optical conductivity has an absorption edge at 
u> . The coherence factors are affected by the broken lat- 
tice inversion symmetry in the CDW phase. The /—sum 
rule reveals an anomalous suppression of the diamagnetic 
spectral weight, mainly for // = 0. This behavior is an 
evidence that there are missing high energy diamagnetic 
states in the SC phase, which would be attracted from 
the bottom to the top of the lower band by the opening 
of the SC gap. Close to the normal-SC phase transition, 
however, these states can by introduced by the same pro- 
cedure that fixes the anomalous divergence of the London 
kernel in the ultraviolet, which is due to the absence of 
the Brillouin zone periodicity into the calculation. We 
have extended the calculation to a general CDW band 
where the loss of the crystal inversion symmetry is in- 
cluded by assumption. 

In summary, we have presented a complete theory for 
s-wave superconductivity in nodal liquids. We have cal- 
culated the thermodynamics, the various response func- 
tions, and transport properties of this system and have 
shown that these quantities deviate strongly from the 
same properties in ordinary BCS superconductors when 
there is particle-hole symmetry. We believe our theory 
can be applied to some TMD, such as 2H-TaSe2 or 2H- 
TaS2 , and our predictions can be checked experimentally. 
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Appendix A: GAP EQUATION 

In this appendix we derive Eq. l|2T]) — ll23l) . Applying 
the variable substitution v — vpk + a^,, the equation 11311 
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can be written into the form 
gvF sr^ f" ,t k 



1 = 
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<T = ±1 



dfc 

^k,afi 



tanh ft 
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In the |/x|/A s <C 1 limit we find: 



1 = 



9 



Aitvavf (P 
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4 In 



cosh (/3a/2) 
cosh (/3A s /2) 



tanh ( -A, 



that is equivalent to Eq. ffil . We notice, however, that 
the above expression remains valid at T c (i.e. for finite 
/j, and A s — > 0) if the strong coupling approximation 
|/x|/Ao M -C 1 is satisfied. 

We define a = 2-Kvpv^l 'g c , Close to T c , taking A s — > 
we obtain, 
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In 



cosh 



/3 c a 



2a 



The critical temperature for 5 > g c is, 
1 



2fc B In 4 



A 2 + fi 2 i n 4 



(A2) 



where A = A S (T = 0,.g,^ = 0) = a(l - g c /g). The 
expression that gives the critical dependence of the gap 
with temperature for |/z|/Ao M <C 1 follows directly from 
the expansion of the gap equation ltl3|l in terms of (3A S . 

To calculate the critical temperature in the weak cou- 
pling regime, we take A s — > in Eq. (| A If) leading to: 



2 *L = -L In 

g p c a 

x tanh 



cosh(/3 c a/2) 
cosh (p c n/2) 

PcH 



2Jm| 
a 



In 



PM 



PAlAI 2 1 a 

dC 



cosh 2 Q 



after integrating the second term of Eq. IIAlj) by parts. 
If Pc\lA <L 4, the integration above can be extended to 
the interval [0, 00], 



25c 
9 



2| M | 



hi 



PA1A1 



In weak coupling (|/i|/Ao M ^> 1) the condition /3 c |/i| ^ 1 
is easily satisfied. The equation above implies that 



r - ZLp"(i-9 C /9)iMr 1 -i 

J-c — , e 
Kbit 



(A3) 



where 111(7) ~ 0.577 is the Euler constant. 



In the weak coupling regime, we can find the gap equa- 
tion near the phase transition. For P\fi\ ^> 1, we use the 
power series expansion in pA s <C 1 of the integral^ 



n tanhff^ 2 + Al 

dv \" 

V^TAf 



dv — tanh ( —v 
o v V2 

7C(3) P 2 Al 



8 7T 2 

Expanding the gap equation l|AlJl in lowest order around 
Pc we find 



A S (T->T C , M ) 



7((3) 



1 



8tt 2 



2&V 



T,-T 



(A4) 

The weak coupling expansions given above are correct 
whenever tanh (fi c \ |/2) ~ 1, or /3 c |/i| > 4. 



Appendix B: SPECIFIC HEAT 

In this section we calculate explicitly the specific heat 
jump in the weak and strong coupling limits. The entropy 
of the problem is given by: 



S = - 



k,7,cr 



-Km ln "k, CT/ x 



where 



^ e 7/3^k. CTM _|_ 2) 1 i s the Fermi-Dirac dis- 
tribution, indexed by a = ±1, and by the two branches 
of the cone. 7 = ±1. The specific heat yields^ 

dS 



C v = -P 



dp 



-k B p 1 



PdAt 



dE k ,^ V 2 dp 

At the phase transition, the specific heat jump reads: 

PI dA 2 



(Bl) 



AC(p c ,n) = lim 

P^Pc 



E 



2 dp ^ ' dE K a^ 



k B Pl dA 2 



8ttvavf dp 



Pc 



E 

<7=±1 



de e sech 



2 / /3c(e + cr/x) 



If P c a > 4, we may extend the integration range to in- 
finity. This integral can be evaluated in two limits, for 
P c \n\ < 1 and p c \n\ > 1: 

k B Pc dA 2 



AC(Am) 



2ttvavf dp 

Pin 



Pc 



ln4+a£- , for/3 c |/*| « 1 
P c \fi\ , for/3 c |/i|»l 



(B2) 
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From Eq. {221 and O, we find: 

^(l n4+ ^) , for/3 c | M |«l 



dA 



d0 



+ -for >' /' I 



(B3) 



In the normal phase, the specific heat CV„ is obtained 
from Eq. jBlt . 



Cvnific) 



xsech 

Evaluating the integral gives: 
k B 1 



pa 

/ dee(e + ^) 2 
„ Jo 



c (e + cr/i) 



18C(3) , \n\0 c « 1 



2nv F v A ft \2 2f3M ;H/3c>>1 



Cy„(/3 C ) 

Combining Eq. llB2l) . llB3|) and l|B4|l . we find 

ACy 



(B4) 



|ln4 A 4+^) , H/? c «l 



2tt 2 tc(3) . i 



H/3 C » 1 



(B5) 



Appendix C: SUSCEPTIBILITIES 

We define the charge and spin susceptibilities from the 
imaginary time ordered correlation functions: 

X c (q,H = -/ dre^ <T T [p(q, r) p(-q, 0)]> (CI) 
Jo 



/ dre^(T r [^(q,r)5 b (-q,0)]), 
Jo 



(C2) 

with p and 5a respectively as the charge and spin density 
operators defined by Eq. l|2"§j) and (|3"4T) . 

The optical, thermal and thermoelectric correlation 
functions are defined as, 

Uij(q,iu)=- dre^OXq,r)ij(-q,0)) (C3) 
Jo 



L ij 



dre^ Of(q,r)jf(-q,0)) (C4) 



iif^ito) =- /Ve^of (q^)ii(-q.o)), 

Jo 



(C5) 



where j is the electric current operator lUfill and j £ is the 
thermal current operator defined by Eq. 1551 . 



Appendix D: HAMILTONIAN IN THE 
BALIAN-WERTHAMER SPACE 

In this appendix we discuss Eq. . The BW space 
is introduced to extend the pairs space (k f, k j) to a 
larger one where the spin and momentum degrees of free- 
dom are decoupled. The procedure rests on "duplicating" 
the Hamiltonian (keeping it invariant by summing in half 
Brillouin zone) , interchange the order of the ip fermionic 
operators in the duplicated term and explore the sym- 
metry under the k — > — k exchange in the k— sum. The 
CDW Hamiltonian in the BW space reads 



Hcdw = J2 ^L(k)k-rf fy 6(T (k) 

V F \ - 



2 ^ 



^ T (k)»pfy 6T (k) 

; (k) rf V, j (k) + ^ |(— k) ff >J r( _ k) 
+V, H (-k)jf a vi i (-k) 
UF^^k) (7aT rf- k*(k) , 



E 

kGiB.Z 



(Dl) 



by the definition of the BW spinor (13311 . 

The chemical potential term © can also be written 
as — ^SkeiBZ ^(k) C0T37/0 *(k). The pairing term 
can also be obtained with the use of the antisymmetric 
property of the Pauli matrix 772 under the transposition 



•" '" " —T]\ a , namely, 



Hp = 



k.a.b 



2 ^ 

k,a,& 



A s V aT (k)% a V b t i (-k) + / l . C . 
^ aT (k)% a V oi (-k) 



A s 



+^a| (-k) % a ViT (k) + 4 1 (k) 4 >1 T (~ k )" 

+^ T (-k) % 6 Val(k)] 

= - £ A s *t (k)(73ri7?2 ^( k) . (D2) 

ke^B.Z. 

Appendix E: LONDON KERNEL 

In this appendix, we evaluate the London kernel 116 71 . 
It can be derived from the calculation of the expectation 
value of the current density operator, 



J(k)=Tr^[e (fiSfCt - f^fl^ck 



(tft( k ) To7?3 tf( k )) 



+e (a,A ck - ^^i^A ck ) (* f (k) r 7?i*(k)) 



(El) 
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in first order in A, where in our definition k = k— f T3A. and 

Expanding the Green function G (iw n ,k) = {iu) n — lo^ 
) _1 up to leading order, 



Tr(*+(k)T^*(k)) = -Tr^r^Go 



1 - - Go 

c 



where Go is the Green function 11211 . The zeroth order 
terms are: 



Tr(* t (k)r r ?3 *(k))c 



7j Tr E ro7 7 3 Go 

x [n(E w )-n(-£J k)(TM )],(E2) 



and 

Tr(* t (k)r r ?1 *(k)>o 



*ck 



/i<7 + El 



E 



x [n(E k<afl )-n(-E^)],(E3) 



where = ^/e£ + A^ k and 

Sk, CTA1 -[(^ + ^) 2 + A2]i 

To first order we find, after a straightforward algebra: 
Tr(*t(k)r ry 3 *(k))i = --A^ - £ Tr [(Go) 2 " 

E 

:1 

[n(-Bk, CT /i) - n(-Ek,ff M )] 



er=±l 

d 



BE, 



a 1 1 



(T=±l 



dE n 



(E5) 



The London kernel IjoTIl follows from the direct substi- 
tution of Eq. (JEU, JUJ, (JEU and (jEHJ into Eq. (JETJ, 
and by noting that the zero order current term, 



U)o = E E e ( e k^i£k + AckgjAck) 
<7=±1 k 

x K-Ek.crp) - n(-Sk j(rM )] 



E£ -Ek.cr/J 



(E4) vanishes by symmetry when integrated over k. 
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